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We study the cosmological evolution of a tachyon scalar field T with a Dirac-Born-Infeld type 
lagrangian and potential V{T) coupled to a canonically normalized scalar field with an arbitrary 
interaction term B{T, <j)) in the presence of a barotropic fluid pt, which can be matter or radiation. 
The force between the barotropic fluid and the scalar flelds is only gravitational. We show that the 
dynamics is completely determine by only three parameters Ai = —Vt/V^^^,X2 ~ —Bt/B^^^ and 
A3 = —B^/B. We determine analytically the conditions for Ai under which the energy density of 

' T, <j> and pb have the same redshift. We study the behavior of T and 4> in the asymptotic limits for 

' A and we show the numerical solution for different interesting cases. 

, The effective equation of state for the tachyon field changes due to the interaction with the scalar 

CN ■ field and we show that it is possible for a tachyon field to redshift as matter in the absence of an 

interaction term B and as radiation when B is turned on. This result solves then the tachyonic 
, matter problem. 

S: 

. I. INTRODUCTION 

. String tachyon fields T are the lowest energy state in unstable Dp-brane or brane-antibrane systems [1]. Since they 

' represent the low energy limit of string-brane models the phenomenology of the tachyon field is important and a great 

. amount of work has been invested in studying the dynamics of tachyon field [1],[2]. In the case of a Dp-brane systems 

1' in string theory, the potential V(T) has been conjectured to be tachyonic at the origin T = [1], i.e. the potential 

O , has a maximum at the origin with a negative mass square, = Vtt < 0. 

, ' The tachyon field has a Dirac-Born-Infeld type lagrangian and therefore it does not have canonical kinetic terms. 

[ This implies that the naive prescription for the mass of tachyon as the = Vtt does not hold in this case [5] nor 

^— I ' can we assume that the evolution of T is to reach the minimum of the potential V{T), as for a standard scalar field, 

K>" ' since V{T) does not correspond to the true potential for T. 

It is well known that the tachyon field has an equation of state parameter — 1 < w < [2]. For string motivated 

potentials, e.g. V ~ the late time behavior gives w = and T behaves as matter. This is the "matter 

^vq problem" because T can easily dominate the universe well before radiation-matter equality since it redshifts slower 
than radiation. So, if T is present at early times it should necessarily decay into other particles to avoid the matter 
problem. 

' In this letter we study the generic cosmological evolution of a scalar field T with a Dirac-Born-Infeld type lagrangian 
, coupled to a canonical scalar field through an arbitrary interaction term i?(T, 0) in the presence a barotropic fluid, 
which can be matter or radiation. We will call the field T as the tachyon field even though we do not constrain 
ourselves to a potential V{T) with < at the origin, so our results are valid for any potential V{T). We show 
that all models dependence is given in terms of three parameters Ai — —Vt/V^^^, A2 — —Bt/B^/^ and A3 — —B^/B. 
^ ' We determine the dynamical equations and obtain the attractor solutions as a function of these Ai parameters. We 
show in this letter that it is possible for a tachyon field to redshift as matter in the absence of an interaction term 
B and to redshift at late times as radiation due to the interaction term. Therefore, the interaction term solves the 
tachyonic matter problem. 

This letter is organized as follows. In section II we set up the framework for the cosmological evolution of two 
scalar fields, a tachyon and a canonical scalar field, with an arbitrary potential in the presence of a barotropic fluid. 
In section III we derive the dynamical first order differential equations and we show that the system is determined 
by only three parameters. In section IV we calculate the critical points. In section V we study different asymptotic 
limits and we present a discussion on specific particle physics motivated examples in section VA. Finally in section 
V B we give some interesting examples and we present our conclusions in section VI. 
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II. COUPLED TACHYON AND SCALAR FIELD 



Our starting point is a universe filled with two scalar fields T, </> and a barotropic energy density ph, which can be 
either matter W{, = or radiation Wb = 1/3. We will assume that the scalar fields interact via a potential B{T,(j)) 
while there is only gravitational interaction between these fields and the barotropic fluid. This work generalizes that 
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of a single scalar field and a barotropic fluid [3], a tachyon field and a barotropic fluid[6] and two single scalar field 
with arbitrary potential and a barotropic fluid [7] . 

One of this scalar fields, namely T , is a tachyon field which is motivated by string theory and D-branes, while the 
other scalar field ^ is a standard canonical field. We take the following Lagrangian for the scalar fields cj) and T [1] 

L = -V{T)^l - d^Tdi^T + ^d^<pd^<P - B{cf>, T) (1) 

where the tachyon is given by a Dirac-Born-Infeld type lagrangian with a potential V{T) and B{(j),T) corresponds 
to the interaction between </> and the tachyon T. The potentials V{T) and B{^, T) are completely arbitrary. The 
equation of motion of (j) and T for a spatially flat Friedman-Robertson- Walker (FRW) universe are 



> + 3H,p + B^ = (2) 



T ■ Vr B-r I — 

+ SHT+-^ + -fVl-T^ = (3) 

where the subindex in V and B is defined as Vr = dV/dT, B^ = dB/d(p and Bt = dB/dT. The Hubble parameter 

H = d/a is 

SH"^ = p = + PT + pb (4) 
where we have taken SttG = 1 and p is the total energy density, pb the barotropic fluid and p,f,,pT are deflned as 

P4> = l^^ + V {<!>), (5) 

for (j) and 

PT = -7^=, PT = VVl - T2. (6) 

\/l-T2 

for T with Pfj„PT the pressure of ^, T, respectively. We deflne the ratio of energy densities as = PtI'^H^-, 
^0 = p^/'iH'^. Using eqs.(5) and (6) we can rewrite the dynamical eqs.(2) and (3) in terms of the energy densities as 

P0 + 3i?p0(l + W0) = BTf = 5 

PT + iHpril + wt) = -Bt f = -6 (7) 
Pb + 3Hpb{l + Wb) = 

where we have included the evolution of the barotropic fluid pb and 

SeeBtT (8) 
defines the interaction term. The equation of state parameters are given by 

Wd, = — = ^. , Wt = — = -l + T . (9) 

* P<p + ' p^ ^' 

In order to have a real energy density for the tachyon we require < < 1 and from eq. (9) we see that the equation 
of state parameter for T is constraint to — 1 < wt < 0. The time derivative of H is given by 

H = -^{p^ + PT + Pb+P^+PT+Pb) = -^(j>'^+ PtT^ + P6(1 + Wb)) ■ (10) 

A. EfFective Equation of State 

From the last two terms in eq. (3) we can define a T-derivative of an effective potential 14// as 
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where we have used px = V/ yl — T^. The dynamics of the tachyon field gives a vanishing dV,,ff/dT, i.e. dVej f /dT 

\ 1 + BtV 1 — )V/p^ = 0. It is important to include in eq.(ll) the multiplicative factor (1 — T^). The dynamics 
of a scalar field is to minimize the potential (i.e. the derivative of the potential w.r.t. the scalar field should vanish) 

and in the absence of the interaction term B the dynamics of the tachyon field does not give a vanishing Vr/V but it 

has a vanishing (1 — T-^)Vt/V. For example for a typical tachyon potential V ~ one has Vr/V = —T and the 

solution to eq.(3) gives t"^ = \,T ^ oo with Vt/V -oo but (1 - f^)VT/V = VrV/p]. = 0. This shows that Vt/V 
is not the T-derivative of the potential for the tachyon T. 

At vanishing dVeff/dT, eq.(3) becomes f + 2,Ht{l - T^) = -dVeff/dT = and it has two solutions: 

V = 0, f = 0, = 1 (12) 

and 

Vt = -St Vl - T2, , ^ = k a{t)-^ (13) 

VI -r2 

with k an integration constant. The solution in cq.(12) necessarily implies V = since nr = y/(3iJVl - 7^2) must 
be smaller than one while eq.(13) gives at late times T = and Vr = —Bt- 

In order to define an effective equation of state it is convenient to rewrite eqs.(7) as 

= -3Hp^{l + W4,eff), PT = -3Hpt{1 + WTeff) (14) 

with the effective equation of state defined by 

BtT BtT 

We see from eqs.(14) that w^eff, WTeff give the complete evolution of p^ and pt- For BtT > we have Wgff < 
and the fluid will dilute slower than without the interaction term (i.e. BtT = 0) while pT will dilute faster since 
WTeff < w^. Which fluid dominates at late time will depend on which effective equation of state is smaller. The 
difference in eqs.(15) is [3] 

AtWe// = WTeff - W4,eff = Aw - T (16) 

with Aw = wt — w^ and T defined as 



= ^ ( ) = ^ ( ) = ^ ( ) (17) 



f P4> + Pt\ 


^^BtT 


( P4> + Pt\ 




/O^ + SItA 


V P't'PT ) 


' m ' 


\ P4>PT J 




[ n^fiT J 



while the sum gives 



flTWTeff + ^(pW^eff = ^TWt + ^4,W^. (18) 

Clearly the relevant quantity to determine the relative growth is given by T and if T < Aw wc have Aweff > 
(i.e. WTeff > ^4>eff) and p^ will dominate the universe at late times with 17^ = l,f^T = 0. For T > Aw we have 
Aweff < and WTeff < w^eff with Pt prevailing and Ot = 1,^^^ = 0. If both fluids have the same redshift, i.e. 
Aweff = and Weff = WTeff = w^eff eq.(18) gives 



Or _ W^, - Weff 
Weff - Wt 



(19) 



and for w^ < wt we have 



while for w^j, > wt we get 



WT^T + w^fl^ 

W4, < W^eff = WTeff = Q ^ '^^^ (20) 



WT^h + w^n^ 

W4, > W^eff = WTeff = + ^ (^l) 

i.e. the effective equation of state is constraint between w^ and wt- In the limit of no interaction 6 = B^^ = we get 
T = and Aw^ff = Aw and depending in which term is smaller w^ or wt we will have either p^ or pT dominating 
the universe at late times. 
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B. Canonical Scalar Field 



In the limit of having only a canonical scalar field (p and a barotriopic fluid, the potential B depends only on (p. If 

the vacmini expectation value of 4>,nin ha-s a finite value, i.e. the potential is of the type B = hgcff^ , then the scalar 
field (j) oscillates around its vacuum expectation value (v.e.v.). If the scalar field has a non zero mass or if the potential 
B admits a Taylor expansion around (pmin then, using the Hopital rule, one has limt_>oo|-B^/-B| = oo and the energy 
density redshifts with = {n — 2)/(n + 2), i.e. w^f, = 0, 1/3 for n = 2, 4 [3]. On the other hand, if (j)min = then 
(j) will not oscillate and \B^/B\ will approach either zero, a finite constant or infinity. Only in the case \B^/B\ going 
to zero or a constant smaller than \/2 will the universe accelerate at late times [3] . 



C. Tachyon Scalar Field 

If we have only a tachyon field in the presence of a barotropic fluid, the late time attractor solution were studied 
in [6]. A tracking solution with constant Ai = —Vr/V^^'^ is given for a potential V = Vo/T^- In this case one finds 
T = yiXi/V^ and yj = V/3H^ = (VAf+36 - A?)/6 with u;t = -1 + = -1 + ( ^/Af + 36 - A?)A?/18 if Wb > wt 
and Wt = wt, if wt, < —1 + (V A| + 36 — Xi)Xl/18 with {xi,yi) = {jb, ±^3^b/Xi). [6]. Clearly a |Ai| <C 1 gives an 
equation of state wt ~ — 1 and an accelerating universe. 



III. GENERIC DYNAMICAL ANALYSIS 



To determine the attractor solutions of the differential equations given in eqs.(2) and (3) or (7) it is useful to make 
the following change of variables 



XI ^ f, y,^l-^Y^ (22) 



and eqs.(7) and (10) become a set of dynamical differential equations of first order 



a;ijv = ~{l-xi)[ 3x1 - V3 Aiyi - V^(l - xl) A2 ^ 

Hn \/3 . 2 ir,A\ 
yiN = — ^ 2/1 - — Ai Xl yi (24) 

(Hn \ / 3 , 2 

3 + -jj- I a;2 + Y - A3 ^2 

Hn ^ \/3 , „ 2 /3 



y2N = Jj- 2/2 - ^ A2 Xl ^2 - Y 2 '^3 X2 y2 

H 3 3 

= -- {Qi^i + + ^b^b) = (xlfli + 2x1 + Qblb) (25) 

£1 Z i 

where N is the logarithm of the scale factor a, N = ln{a), 76 = 1 + Wb, 71 = 1 + wi, 72 = 1 + W2 and /at = df/dN 
for / = H,Xi,yi {i = 1,2), Qb = ^ — ^1 — ^2 and 

^i(^)--^' ^2(iV)--;^' HN)^-^. (26) 

Notice that all model dependence in eqs.(24) is through the three quantities Xi{N),i = 1,2,3 and the constant 
parameter 7^ = 1 + Wb- The last equation of (24) is constraint between — 3 < Hn/H < for all values of Xj, and 7f,, 
it takes the value —3 when the universe is dominated by the kinetic energy X2 = 1 and therefore fib = fli = while 
it becomes Hn / H = when X2 = Xi = Oft = and the universe is dominate by a constant potential 2/1 + J/| = 1 . 
The set of equations given in eqs.(24) give the evolution of two scalar fields <j)^T^ a tachyon and a canonical scalar 
field, with arbitrary potentials in the presence of a barotropic (perfect) fiuid with equation of state = 1 — 7b- If 
we do not want to consider the contribution from the barotropic fluid we can easily take the limit 7b = in eqs.(25) 
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since all contribution form is given in Hn /H via the term f7b7(,. For fi^, 7^ we will assume a barotropic fluid with 
< 76 < 2 and 76 = 1 for matter while 76 =4/3 for radiation. 
In terms of a;,, y, we have 

fi* = n2 = ^ = x^ + !/l, ^^=xl-yl (28) 



and 



Pt , , 2 _ P<j> xl-yl 



W\=Wt = — = -1 + Xi, = Wa> = — = -1) n (29) 

Pt P4> a + 2/2 

with 71 = 1 + wi = xf , 72 = 1 + u'2 = 2x|/ri2- The interaction term defined in eq.(8) and (17) are now 

b = BtT = -■i^l'^H^ \2x1yl (30) 
Y ^ A2Xi?y| / Oi + »2 \ ^2^^ 



V3 V ^2 

giving an effective equation of state parameters defined in eqs.(15) as 

A2a;iy2 \/lVliwi- \2x1yl 
wuff ^ wreff = - = -^^ (32) 

, A2 V3il2W2 + A2 xi y| 

W2.ff ^ W,.ff =W2+ = (33) 

and 

71 eff = 1 + Wieff, 72 e// = 1 + W2ef f , Ibeff = lb- (34) 

The acceleration of the universe is given by 

- = H^{1 + ^) = (^6(1 + 3^6) + Oi(3a;? - 2) + Axl - 2yl) (35) 



where we have used H = HH]^ and cq.(25). Clearly acceleration will occur if the universe is dominated by the 
potential y| = B/3H'^ or Oi with < 2/3, i.e. for yj + rii(l - 3x^/2) > il6(l + 3^6) + 4x|. 

IV. CRITICAL SOLUTIONS 

To find the critical solutions to the dynamical equations (24) we need to solve them for xijv = X2Ar = yiN = 2/2JV = 
for constant values of Aj, i = 1, 2, 3. However, the set of equations are highly non linear and it is quite complicated to 
find all critical points. Instead, we will calculate the critical points when all three fluids have the same redshift. We 
will also study different limits and we solve eqs.(24) numerically to show the behavior of the scalar fields T and (p. 

As discussed in section II A the energy density pi,p2, Pb with the smallest effective equation of state, or alternatively 
with the smallest eff = ^ + Wi^ff, i = 1, 2, 6, will dominate the universe. For ji eff <Jkeff, Ij eff we have ili — 1 
and rife = fij = at late times. If two or the three energy densities have the same (smallest) effective equation of 
state then two fi's or all three Vt's will be different than zero. Therefore, eq.(25) gives asymptotically 

2_/?7v 

= (^l7l + ^^272 + ^blb) = (^l7l eff + ^2l2 eff + ^blb) = lef f (36) 

where the last equality in cq.(36) is valid only asymptotically and 7e// = 1 + w^ff is the smallest of the "y^ff^s, i.e. 
7e// = smallest{ji eff,l2 eff,lb)- 

If one of the scalar fields dominates over the other one then eqs.(24) will be reduced to a single scalar field, either 
T or (/), in the presence of a barotropic fluid. These special case have been studied in the literature for tachyon 
domination [6] and for canonical scalar field in [3]. 
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FIG. 1: We show for Ai — 10, A2 = —10, A3 = 5 and 'yt = 1 + Wb = 1 the evolution of Qi = ^t,^2 ~ Q,p,ilb (blue (solid), 
red (dotted) and black (dashed), respectively). We also show the equation of state parameters wi = WTjWTeff (blue, yellow, 
respectively) and W2 = w^jW^^ff (green, red, respectively) as a function of A'^ = Log[a]. With these choice of A's the attractor 
solution has {xi,yi) = (0.48,0.36) and (0:2, 3/2) = (0.48,0.34), = 0.14,^2 = 0.35,^6 = 0.51 and wieff = W2eff ^ wt = 0. 
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FIG. 2: We show for Ai = 50, A2 = 0, A3 = 3 and 7;, = 1 + u;;, = 4/3 the evolution of 57i = Q.t,^2 = Oi^,f2f, (blue (solid), 
red (dotted) and black (dashed), respectively). We also show the equation of state parameters wx = wt = Wxeff (blue) 
and W2 — — w^^jj (red (dotted)) as a function of = Log[a\. With these choice of A's the attractor solution has 
(a;i,yi) = (0.99,0.03) and (3:2, 2/2) = (0.40,0.40), = 2/3,^2 = 1/3,^6 = and wie// = W2eSf = 0. 



So, let us now determine under which conditions do both fields T and 4> have the same redshift. We see from 
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FIG. 3: We show for Ai — 50, A2 = —3, A3 = 3 and 75 = 1 + = 4/3 the evolution of f2i — Qt, ^^2 ~ ^(/i, (blue (solid), 
red (dotted) and black (dashed), respectively). We also show the equation of state parameters wi — WTjWTeff (blue, yellow, 
respectively) and W2 = w^,w^^]j (green, red, respectively) as a function oi N = Log\a\. With these choice of A's the attractor 
solution has (xi,t/i) = (0.28,0.16) and {xi,yi) = (0.63,0.41), fii = 0.03,^2 = 0.56,^6 = 0.41 and wxeff = W2e// = wi, = 1/3. 
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FIG. 4: We show for Ai — 100, A2 — 1/2, A3 = 1 and 7;, = 1 + u;;, = 1 the evolution of Q.\ = 57T,f22 ~ ^^■,Q.b (blue (solid), 
red (dotted) and black (dashed), respectively). We also show the equation of state parameters lUi = WTjWTeff (blue (dotted), 
blue (solid), respectively) and W2 = w^yW^efj (red (dotted)), red (solid ), respectively). With these choice of A's the attractor 
solution has {xi, j/i) = (1, 0) and (2:2, j/2) = (1/3, 0.78), fii = 0.27, ^2 = 0.73, Q.b — Q and wieff = W2eff — Wb = —0.5. 



eqs.(32) that both effective equations of state are equal if 



(37) 
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FIG. 5: We show for Ai = 1, A2 = 1,A3 = —1 and 71, = 1 + uij = 1 the evolution of f2i = Q.t,^2 = f20,r26 (blue (solid), 
red (dotted) and black (dashed), respectively). We also show the equation of state parameters w\ — WT,WTeff (blue, yellow, 
respectively) and 1112 = w^jW^^ff (green, red, respectively) as a function of = Log[a\. With these choice of A's the attractor 
solution has (2:1,1/1) = (0.55,0.8) and (x2,t/2) = (-0.11,0.49), fii = 0.74,^2 = 0.26,^6 = and wie// = W2e// = -0.75. 



and using eq.(36) and eq.(37) we get 



^1 _ 71 - 7e// _ Wl - We// 
^2 leff - 72 Weff - W2 



(38) 



which coincides with eq.(19). The solution of the dynamical equations given in eqs.(24) for finite values of X2,yi, 2/2 
as a function of a;i are 



yi 



X2 



2/2 



Xi 



\/37, 



Xixi 
2" Aayi 



3 (2 - 7e//) 
/ 3V37g^^(x^-7e//) 

1^ xlx^xl^T^, J 

2 /37e// A2y2 



1/3 



G(a:i) 



V3A22/2 V 2 2-7e//^ 
for yiAT = 0, X2N — 0, xiat — 0, j/27V = 0, respectively. We have defined 



1/3 



and we obtain 



37f 



xlxj^/T^, 



(39) 
(40) 

(41) 
(42) 



(43) 



(44) 
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where we have used cq.(39) in eq.(44) and (41). The values of 2/1,0:2,2/2 are functions of xi and 7e//- Sohitions with 
fii ^ 0, J72 7^ 0, fit, 7^ imply that "feff = 7fc — li eff = I2 eff is a given constant (given by the barotropic fluid). 
If however, 7e// = 71 eff — 72 eff < 76 then fJ;, = and 7e// is not a constant given a priori but it is a constant 
which is a function of Xi,yi. In this case eqs.(39)-(42) must be supplemented with eq.(36). Using eqs.(39)-(41) we 
can rewrite eq.(36) as 



3. -12 5/3 

. + \. '^iJ G{x,f (45) 



^i7i _|_ ^272 _ 37e// ^ 3 4A37g^ ^ 
7e// leff ~ Xl^/l -xl (2-7)^ 



Ai^^A2^'^(2~7e//)a ^ ,2 . ^leff{'i-leff)a r.QS 

^ 1/3 > ^3 > 5 (49) 



which is only valid if fl;, = 0. 

Let us take 7e// = 76 with a barotropic fluid with 76 > 1, i.e. it includes matter 76 = 1 and radiation 7b = 4/3, 
and in this case we get the following constraints on Ai, A2, A3, from requiring |a;j| < 1 and < yi < 1: 

|A?A2| > |^^!2|ij!£_fil| > 3V37.^,,Gf„ (47) 
Mxi < 0, Xixi > 0, A3a;2 > 0, A1A2 < (48) 

^ r 

^llif-Gl 

where Gm in cqs.(47) and (49) is the minimum value of G{xi). The constraint (46) comes from requiring fii < 1 
and (xf a/1 — xf)^^ has a minimum value (xj^l-xj)-^ = 3^3/2 at xj = 2/3. Eq.(47) comes from eq.(41) and 
< 2/2 < 1- Eq.(48) arises from requiring < yi and < 2/2 and we see from eq.(41) that GA2 must be negative, 
which for 7e// > 1 implies that A2.T1 < 0. Finally cq.(49) comes from eq.(50) and < 2/2 < 1- 

For 76 = 7e// = 1 the function G is simply G = (1 — x^Y^^/xi and G takes the value |G(0)| = 00 and G(±l) = 
Gm = so that < |G(.x)| < 00. Therefore, the lower and upper constraints in eq.(47) and (49), respectively, do not 
impose any limitation to A^. On the other hand if 7b = ^^ff > 1 (as for radiation 75 = 4/3) then |G(±1)| = |G(0)| = 00 
and G(x) has a minimum absolute value Gm = \G{x — Xm)\ at x"^ = [1 + ^leff — \/ ^ + 167e//(7e// — l)]/2, e.g. for 
7e// = 4/3 we have x„i — 0.93 and Gm = 1.16. The fact that G has a minimum value for 7,.// = 76 > 1 implies a 
constraint on the value of Ai,A2 namely AfA2 > 3\/37gjjG^ as seen from eq.(47) and from eq.(49) an upper value 
for A3 . 

To finally solve eqs.(39)-(42) we substitute 2/2 from eq.(41) into eq.(42) and we get a nonlinear equation for xi only 

4/3 

„2 _ ^"feff ^2 _ 3(2 -7e//) / AaXiyaA _ 3(2 - 7e//)7e//Q: , . 



with 



a=l 



M \ -1/3 



] j;;i'x,G{xx) (51) 



and a depends on A2/A1 and xi with a > 1 since the second term in eq.(51) is positive for 7e// > 1- Eq.(50) is 
however valid for any value of 7(,j/ including 7ej/ < 1. In this later case the second term in eq.(51) is no longer 
positive definite {xiG is negative for > Jeff) and the fact that a must be non negative since 2/2 > 0, gives the 

constraint X2^^ Xi^^'^j^^j^xiG{xi) > —1. 

The solution to eq.(50) is non analj^tical and must be solved numerically for xi. In terms of G eq.(50) gives 

, 4/3 > 2/3 N 

G{xxr = ^^^^^^>Gl. (52) 

For 7b = 7e// = 1 eq.(52) does not impose a constraint on the choice of Ai, A2, A3 since G can take values form zero 
to infinity and Gm = 0. However, for jeff > 1 the function G has a minimum value given by Gm, therefore the choice 
of Ai, A2, A3 must be such as to give a value of G larger than Gm and for jh = 4/3, as mentioned before, we have 

Gm = 1.16. 
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Regarding the quantity a given in eq.(51) which term dominates (first or second) depends on the vahie of A2/A1 
and xi through the function F{xi) = xiG{xi) = {{'^eff — — xf y^^. The function F has a minimum value at 

Xi = 2 — Jeff giving a value Fj^ = 2^/^(7e// — 1)^/^. For 7e// = 1 one has < F = ^/l — x\ < 1 while for 7e/j^ > 1 
one finds < F < F{x = 1) = 00. However, since F(0) = F{xi) = 7^^^ for xi = \/jeff^ — leff) the value of F 
in the case 7e// > 1 is constraint between Fm < F < F(0) = F{xi) = 7^^^ for < \x\\ < ^7e//(2 — Jeff), e.g. for 
Jeff = 4/3 we have F„ = 1.05 < F < F{0) = {4/3)^/^ = 1.1 for < |a;i| < 2^2/3 = .94 which covers most of the 

— 1/3 

range |xi| = (0, 1). Therefore, Tg// ^ is in general of order unity giving a in eq.(51) as 

a ~ 1 + (A2/Ai)2/^ (53) 

If we have jeff = 7i e// = 72 e// < 7b then fif, = and jeff is not a given constant but is a function of Xi, jji and 
may be smaller than 1. In this case the solution to eqs.(39)-(42) is given by solving eqs.(45) and (52) numerically. 
If Jeff < 1 then the minimum of \G\ and |F| are Gm = F^ = with < G < 00 and < F < 00 and there is no 
upper constraint on A3 from eq.(49) and no lower constrain on IAJA2I in eq.(47). However, since F(xi) = XiG{xi) = 
{{jeff — — x\ can now be negative for x\ > jeff the condition 2/2 > (i.e. a > 0) implies from eq.(51) 

the constraint 



leff - xl Xjjef f 



F{x,r = ^l=^ > (54) 



V. ASYMPTOTIC BEHAVIOR 



Special cases can be studied by taking different limits of the parameters Ai. A constant Ai is given by a potential 
V = Vo/T'^ and wc have Ai = 2/v^l4- A constant A3 can be obtained for example if we consider a factorizcable 
interaction potential B{(j),(f) = h{T)K[(l)). In this case A3 = —K^/K and taking K{(l)) as an exponential potential 
we have A3 = — /3 for K = KqC^'^. However, a factorizcable potential B gives A2 = — /it/(/i'^^^^^''^) which is now a 
function of both fields T and So, a constant A2 can only be obtained as a combination of the evolution of T, <j) and 
cannot be trivially anticipated just by looking at the functional form of B. 

An interesting limit is when x\ = T'^ = 1. In this case j/i = since Oi = yl/ ^Jl — x\ must be smaller than one. 
The value of Vli depends on the value of the parameters Aj. Solving eqs.(24) we find two solutions. One solution has 
x\ = l,yi = and 

X2 = (^(2 - 7e//)Ai + 47e//A^ - \2sgn[xl]^J^,{2 - Jeffheff^l + (2 - lefffxi^ 

" ^§ (~^^ ~ 7e//)A2S5n[a;i] + ^8(2 - jeff)leff\l + (2 - lefffxi^ (55) 

where sgn[x-i\ = ±1 and eqs.(55) constrains the values of A2,A3 since \x2\ < 1,0 < t/2 < 1- We see that this solution 
requires IA2/A3I < 1. The second solution has x^ = l,yi = and {x2,y2) = (0,0), i.e. CI2 = 0, which reduce to the 
case without canonical scalar field discussed below. 

In the limit A2 = the eqs.(24) reduce to two uncoupled scalar fields with a vanishing interaction term 5 and 
71 e = 71 , 72 e// = 72 • The solution can be obtain from the combination of the work given in [6] and [3] , where a 
tachyon field and a scalar field in the presence of a barotropic fluid, respectively, are studied. Which term dominates 
at late times depends on which equation of state is smaller a nd this depends on the values of Ai,A3. The equation 
of state for the field T is given by 71 eff = 7i = (\/Af + 36 — Af )Af /18 if jb > ji and 71 eff = 7i = 76 otherwise 
[6] while for the field (f) we have 72 e// = 72 = A§/3 for A3 < 3jb and 72 e// = 72 = Jb if A3 > 37f, [3]. Therefore if 
A§ < (VAt + 36- A?)A?/6 then we have 72 < 71 and the energy density will dominate over pT at late times. For 
A| > (-^/Af + 36 — A^)A^/6 we have 72 > 71 and the energy density px will dominate over p^. 

In the limit A3 = we have from eqs.(24) X2n/x2 < which implies a decreasing X2 with an asymptotic limit X2 = 0. 

In this case eqs.(24) or equivalently eq.(50) the solution requires a = 0, i.e. F{xi) = xiG{xi) = — (Ai/A2)^/^7gj^. 

Once Xl has been obtained from a = wc can use cqs.(39) and (41) to obtain the values of yi and 2/2- 

Finally, if we take Ai = then from j/ijv = in eqs.(24) we have yiH^/H = 0. If Hn/H = then X2= x\=Q.b = Q 
and from X2n = we conclude that 2/2=0 and the universe is dominate by a constant tachyonic potential Q.i= y1 = 1. 
If however, j/i = then from xin = wc need x\ ~ \ or y2 = 0. For x\ 7^ 1 the solution has j/i = j/2 = 2:2 = with 
Vti = Q.2 = Q and fif, = 1. li x\ = 1 the solution is yi = and X2, 2/2 are given in eqs.(55). 
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A. Particle Physics Model 

Let us now take a specific choice of potentials motivated by particle physics. We consider a factorizeable interaction 
potential 

B{^,^) = h{T)K{^). (56) 
In this case the Aj parameters become only functions of a single field 

^^-~V3/^' ^2-- ^3/2^1/2 ' ^^'^^ 

We see from eq.(57) that Ai is only a function of T, A3 is only a function of <p while A2 depends on T and </>. Let us 
take from eq.(ll) dVeff/dT = {Vr + Bt Vl — )V/ — 0. This equation has two solutions given in eqs.(12) and 
(13). If T 7^ 1 then using the anstaz in eq.(56) we have from eq.(ll) 

m)- — ^ = (58) 

where we have used px = V/ \J\ — T'^ and substituting K in eq.(57) we get 



^^=lt) K-vi) '^«> 

The solution in eq.(13) has T ~ and px — V. As an example let us take the potential for brane-tachyon field as 
V = Vo e-'^'/2 and the interaction term as B = Boe<^'^(l)"', i.e. h = e'^'^ and K = Bocj)'^, with Vo >0,Bo> 0. In this 
case we have from eq.(58) K = pTTe~^'^ /jS 

and in the limit T — > 00 we have 

Ai -^00, A2 ^ —00, A3 — > ±00 (61) 

where we have taken /3 > as required by eq.(ll) since Vr and Bt must have opposite signs. The relative growth is 
given by 

^ ^ 2^3/2 ^ ^ y-(n+l)/rig/3T+T=(2-n)/4n ^ ^ y(n-2)/ng/3T+T=(2-n)/4n j-g2-) 

A2 Ai A2 

Form eq.(62) we have Ai » [A2I and for n = 2 we get IA3I » Ai » IA2I while for n > 2 we find Ai » IA2I ^ jAaj. On 
the other hand if = 1 then wc have V — and the solution has either = or it is given by eq.(55). Solving 
numerically eqs.(2) and (3) with n = 4, /? = 1, 7{, = 1 we find t"^ ^ I with |Ai| > IA2I > IA3I > 1 and ^2 = il.b = 0. 



B. Examples 

We now present four different attractor solution depending on the values of Ai, A2, A3. We show in figures 1-5 the 
evolution of = rii,r2i^ = ^2,^b and wi = w^,W(i)eff,w2 = w^,w^eff for the different choices of X's. We also 
show the phase space of {xi,yi) and (x2,y2) for each case. Since the phase space depends on four variables, namely 
(xi, 2/1,0:2,2/2) , it is no surprising that the curves in the two dimensional space {xi,yi) and {x2,y2) may cross. 

In figs.l we have Ai = 10, A2 = —10, A3 = —5 and 76 = 1 + W6 = 1. The conditions in eqs.(46)-(49) are therefore 
satisfied and we expect to have Vli,Vl2T^b different than zero. In fact wc get Vli = 0.14, J72 — 0.35, fib = 0.51 and 
{xi,yi) = (0.48,0.36) and (^2,2/2) = (0.48,0.34). The effective equations of state are Wieff = W2eff = Wb = giving 
a decelerating universe. 

In figs. 2 we show a model with vanishing interaction term (5 = 0, i.e. A2 = 0, and with Ai = 50, A3 = 3 and 7^ = 
l + Wb = 4/3. The attractor solution has fii = 2/3,^2 = 1/3,^6 = and ixi,yi) = (0.99,0.03) (^2,2/2) = (0.40,0.40). 
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The equation of state Wi is equal to the effective equation of state Wieff = Wi = and T redshifts as matter at late 
times. 

In figs. 3 we show the same model as in figs. 2 but with an interaction term S ^ 0. Wc take A2 = —3 and Ai = 50, A3 = 
3 and 7fc = 1 + lOf, = 4/3. The attractor solution has fii = 0.03,1^2 = 0.56,17b = 0.41 and = (0.28,0.16) 

{x2,y2) = (0.63,0.41). The effective equation of state is given by wie// = W2eff = Wb = 1/3 which shows that the 
tachyon field and <f> redshift as radiation. We see that the relevance of the interaction term in figs. 3 which makes 
Wieff go from Wieff = (no interaction) to Wieff = 1/3 (with interaction). This shows that an interaction term can 
solve the problem of the tachyon field redshifting as matter and dominating the universe well before radiation-matter 
equality. 

In figs. 4 we have Ai = 100, A2 = 1/2, A3 = 1 and jb = 1 + Wb = I. With these choice of A's the attractor solution 
has {xi,yi) = (1,0) and {x2,y2) = (1/3,0.78), Oi = 0.27,^2 = 0.73, = and wie// = 1^26// = Wb = —0.5. Notice 
that even tough yi ^ wc have xi ^ 1 and a finite f2i which dominates the universe with wie// = W2eff = —0.5 
giving a positive acceleration. Notice that ^ as for a tachyon with no interaction, however the interaction term 
gives a negative effective equation of state. 

In figs. 5 we have Ai = 1, A2 = 1, A3 = —1 and jb = ^ + Wb = I. The conditions in eqs.(46)-(49) are therefore not 
satisfied and we do not expect to have Oi, f22, fib different than zero. In fact we get = 0.74, ^2 = 0.26, fib = 0, i.e. 
the T field dominates at late times, and {xi,yi) = (0.55,0.8) (0:2,2/2) = (—0.11,0.49). The effective equations of state 
are wie// = W2eff = —0.75 giving an accelerating universe. 

VI. SUMMARY AND CONCLUSIONS 

We have studied the dynamical system of two scalar fields, a tachyon and a canonically normalized field, with 
arbitrary potentials in the presence of a barotropic fluid in a FRW metric. We have shown that all model dependence 
is given in terms of three parameters, namely Ai = —Vt/V^/^,X2 = —Bt/B^/^ and A3 = —B^/B. The solution to 
the dynamical equations given in (24) are non linear and general analytic solution does not exist. If one of the scalar 
fields dominates over the other one then eqs.(24) will be reduced to a single scalar field, either T or 0, in the presence 
of a barotropic fluid. These special cases have been studied in the literature, for a tachyon domination in [6] and for 
a canonical scalar fleld in [3]. We have determined the critical points in the case where the two scalar fields T and 
(j) have the same redshift. The solution is a function of the parameters Ai and we have calculated the restrictions on 
these parameters. We have solved numerically for different interesting cases and we have shown that either scalar 
field or the barotropic fluid can dominate at late times. The interaction term between the tachyon and the scalar 
field changes the effective equation of state for T and it is possible for a tachyon field to redshift as matter in the 
absence of an interaction term B and as radiation when B is turned on. This results solves then the tachyonic matter 
problem. 
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